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ABSTRACT 

We present high-resolution optical velocity fields from DensePak integral field spectroscopy, along 
with derived rotation curves, for a sample of low surface brightness galaxies. In the limit of no baryons, 
we fit the NFW and pseudoisothermal halo models to the data and find the rotation curve shapes 
and halo central densities to be better described by the isothermal halo. For those galaxies with 
photometry, we present halo fits for three assumptions of the stellar mass-to- light ratio, T*. We find 
that the velocity contribution from the baryons is significant enough in the maximum disk case that 
maximum disk and the NFW halo are mutually exclusive. We find a substantial cusp mass excess at 
the centers of the galaxies, with at least two times more mass expected in the cuspy CDM halo than 
is allowed by the data. We also find that to reconcile the data with ACDM, ~20 km s _1 noncircular 
motions are needed and/or the power spectrum has a lower amplitude on the scales we probe. 
Subject headings: dark matter — galaxies: fundamental parameters — galaxies: kinematics and dy- 
namics 
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1. INTRODUCTION 

The behavior of cold dark matter (CDM) on galaxy 
scales has been discussed extensively in the litera- 
ture, with particular attention given to constraints 
from rotation curves of low surface brightness (LSB) 
galaxies. Numerical simulations of CDM predict 
cuspy halos with a density distribution showing 
a p ~ r~ a (a > 1 ) beh a vior, regardless of mass 
(e.g., ICole fc Lacevl Il996t iNavarro. Frenk. & Whitel 
19961 119971: lAvila- Reese. Firmani. fc Hernandez! 119981 : 
Klypin et all 120011: iDiemand et all I2005D . However, 
numerous studies of the rotation curves of LSB and 
dwarf galaxies have found the data to be inconsistent 
with a cuspy halo, and instead to be more consis- 
tent with a halo having a nearly constant density 
core: p ~ r~ a (a ~ 0) (e.g.. iFlores fc Primackl 
1991 Ide Blok. McGaugh fc van der Hulstl IIqM 
McGaugh. Rubin, fc de Blokl 120011: iMarehesini et al l 
gOOgl: Ide Blok Bosma. fc McGaugh l l200d iGen tile et al. 
2004 120051: Islmon et al.l 120051 : iKuzio de Narav et all 
20061 : ISpano et al.l 120071 ). This interpretation of the 
observations has been met with skepticism, though, 
as cored halos lack theoretical and cosmological moti- 
vation. Systematic effects in the data (most notably 
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beam smearing, slit misplacement, and noncircular 
motions) have also been u sed to argue against the 
presence of co r ed halo s (e.g.. Ivan den Bosch et al.ll2000l 
Swaters et all l2003al: ISimon et all 120031 : iRhee et all 



2004) . although some authors find the magnitude of 



these effec ts to be insufficient for masking the presence 
of a cusp (|de Blok et al.ll2003f K 

To address these concerns, the most recent obser- 
vations have been made u s ing integral field spectro- 
graph s (e.g., IChemin et all 12004 IGentile et all 12004 



20051: IKuzio de Narav et alJ 120061 : ISimon et all 120051 : 
Swaters et al.l l2003bf ). By their nature, these high- 



resolution two-dimensional velocity fields eliminate con- 
cerns about long-slit placement, and highlight the pres- 
ence of noncircular motions. They also typically probe 
in detail the innermost regions of the galaxies where the 
cusp-core conflict is most severe. Even with these im- 
proved observations, the data remain more consistent 
with cored halos and tend to have concentrations too 
low for ACDM (Kuzio de Naray et al. 2006; Gentile et 

al. 2007; but see Swat ers et al. 2003b). 

In a previous paper (jKuzio de Narav et~aDl2006l here- 
after K06) , we presented high-resolution optical veloc- 
ity fields and derived rotation curves for 11 LSB galax- 
ies, along with pseudoisothermal (core) and NFW (cusp) 
halo fits in the limit of no baryons. LSB galaxies are 
thought to be dark matter-dominated down to small radii 
(de Blok & McGaugh 1996, 1997; Boriello & Salucci 
2001; but see Fuchs 2003) with the light simply pro- 
viding a tracer for the dark matter, so neglecting the 
baryons is not entirely unreasonable. The stellar mass 
contribution in these systems is low, which reduces er- 
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rors involving the uncertainty in the stellar mass-to-light 
ratio, T», and in turn, the isolation of the dark matter 
component. While halo fits are often made under the 
assumption that all observed rotation is due to the dark 
matter, it is not strictly true. In this paper, we present 
ha lo fit s for four assumptions about Y» for the galaxies 
in |K06| with available photometry. We also present ob- 
served DensePak velocity fields, rotation curves and halo 
fits for six new galaxies, along with additional DensePak 
observations, up dated rotation curves and halo fits for 
three galaxies in lK06l 

In § 2 we describe our sample of galaxies and the 
DensePak observations; the data reduction is discussed 
in § 3. The observed velocity fields and derived rotation 
curves are presented in § 4 and the zero disk halo fits 
are presented in § 5. Mass models for the galaxies with 
photometry are presented in § 6 and § 7. A discussion of 
the mass model results is given in § 8. Our conclusions 
are stated in § 9. 

2. SAMPLE AND OBSERVATIONS 

During the nights of 2006 March 28-29, April 6-8, Au- 
gust 28-30 and September 25-27 we observed 17 LSB 
galaxies using the DensePak Integrated Field Unit (IFU) 
on the 3.5 m WIYN 3 telescope at the Kitt Peak National 
Observatory (KPNO). UGC 4325, DDO 64, and F583-1 
were observed during March and April. These observa- 
ti ons au gment previous DensePak observations published 
in IK06I . The 14 targets that were observed in August 
and September were selected pri marily from the Nearby 
Galaxies Catalogue (|Tullvl [i~988) . Selection criteria for 
these galaxies included positions satisfying 18 h < a < 
08 h and 10° < 5 < 50°, inclinations between 30° and 
85°, heliocentric velocities < 3000 km s _1 , and an esti- 
mated Vfiat (approximated by Vfiat ~ 0.5W2o(sin i) _1 ) 
between roughly 50 km s _1 and 100 km s . Addition- 
ally, we selected those galaxies which appear to have dif- 
fuse Ha emission and lack indicators of significant non- 
circular motions (e.g., strong bars or gross asymmetries). 
We also targeted galaxies with previous long-slit rotation 
curves or Ha imaging. 

Because of inclement weather and telescope scheduling, 
we did not have Ha imaging or long-slit rotation curves 
for most of the galaxies in the sample prior to making 
the DensePak observations. Without these data, there 
is no way of knowing how much Ha emission will be 
detected by the IFU until the observations are actually 
made. Thus, our observed sample of galaxies is signal- 
limited; the IFU fibers detected sufficient Ha emission to 
create a useable velocity field in just under 50% of our 
target galaxies. 

The obser ving setup and procedure were identical to 
that used in IK06I . The IFU orientation on the sky and 
the total number of pointings per galaxy were tailored to 
each galaxy so that the critical central regions were cov- 
ered by the DensePak fibers. Each exposure was 1800 s, 
and two exposures were taken at each pointing. To pro- 
vide wavelength calibration, a CuAr lamp was observed 
before and after each pointing. We used the 860 line 

3 Based on observations obtained at the WIYN Observatory. 
The WIYN Observatory is a joint facility of the University of 
Wisconsin-Madison, Indiana University, Yale University, and the 
National Optical Astronomy Observatory. 



nun" grating in second order, centered near Ha, giving 
a 58 km s _1 velocity resolution. The DensePak fibers 
are 3", however, we achieved ~2" resolution by shifting 
the DensePak array by small amounts so that the spaces 
between the fibers were observed. The distances to the 
galaxies in the sample are such that a 3" fiber provides 
subkiloparsec resolution. 

3. DATA REDUCTION 

The dat a we re reduced following the procedure de- 
scribed in K06. Briefly, the observations were reduced 
in IR AF 4 using the HYDRA pa ckage. Night-sky emission 
lines (|Osterbrock et al.lll996h were used for wavelength 
calibration, and velocities were measured by fitting Gaus- 
sians to both the sky lines and the four galactic emis- 
sion lines of interest: Ha, [Nn]A6584, [Sn]A6717 and 
[Sn]A6731. The average error on individual emission- 
line velocities due to centroiding accuracies was roughly 
1.5 km s _1 . The velocity assigned to each fiber was the 
arithmetic mean of the measured emission-line velocities 
in the fiber. The error on the fiber velocity was the maxi- 
mum difference between the measured velocities and the 
mean. Many of these errors were less than 5 km s _1 , 
although a few were as high as ~20 km s^ 1 . If only Ha 
was observed in a fiber, the observed Ha velocity was 
taken as the fiber velocity. Without other lines to de- 
termine the maximum difference, we adopt 10 km s _1 as 
a conservative err or estimate based on experience with 
long-slit data (e.g. iMcGaugh et al.ll200l . 

Using the input shifts at the telescope, individual 
DensePak pointings were combined to construct the ob- 
served velocity field. An Ha flux image of the galaxy 
constructed from the DensePak observations was com- 
pared to an actual Ha image of the galaxy obtained at 
the KPNO 2 m telescope to confirm the accuracy of the 
offsets. The fiber maps were registered to the galaxy im- 
age by the correspondence of flux maxima through the 
fibers to observed features like individual Hn regions. We 
find the accuracy of the fib er po sitions (~0".6) to be con- 
sistent with our results inlK06l. Both large (l"-2") and 
small (~0".7) shifts can be made confidently. 

Rotation curves were derived from the observed ve- 
locity fields usin g the tilted-ring fitting program RDTCUR 
(Bcgcman 1989) following the procedure outlined in lK06l . 
To construct a rotation curve, RDTCUR requires the sys- 
temic velocity, inclination, kinematic center, and position 
angle of the galaxy in addition to the observed fiber ve- 
locities. For a more complete descri ption of R0TCUR as 
applied here, the reader is referred to lK06l 

Because the DensePak data cover the centers of 
the galaxies and probe the regime of solid-body ro- 
tation, neither the galaxy center nor the inclination 
could be determined from the observations. The ve- 
locity field centers were therefore fixed to the opti- 
cal centers of the galaxies, as determined by the cen- 
troid of ellipses fit in the surface photometry, and 
the inclination was fixed to published values 
1988; dc Blo k fc McOaughl HoOft IMcGaugh et al 



Won 



de Blok et al J l200lHde Blok fc Bosmdl2002t Uames et al.1 
2004|) . The systemic velocities were determined by 



4 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA), Inc., under agreement with the 
National Science Foundation. 
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TABLE l 

Properties of Observed Galaxies 





R.A. 


Decl. 


lining 






v hel 


1 wnax 


V 

* max 


P.A. 








C 12000 01 


C T2000 0) 


(m£ig sxcsgc ^) 






(km s^ 1 ) 




(km s ) 




( km s ") 


Refer oi ict 's 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


(12) 


UGC 4325 


08 19 20.5 


+50 00 35 


22.5" 


10.1 


41 


514 


2.9 


110 


52 


9.0 


2,2,2 


DDO 64 


09 50 22.4 


+31 29 16 




6.1 


60 


517 


1.9 


60 


97 


7.9 


2,2,2 


F583-1 


15 57 27.5 


+20 39 58 


24.1 


32 


63 


2256 


5.8 


72 


355 


8.7 


4,4,4 


NGC 7137 6 


21 48 13.0 


+22 09 34 


20.7 C 


22.5 


38 


1669 


3.2 


62 


44 


7.1 


1,10,12 


UGC 11820 


21 49 28.4 


+14 13 52 


23.7 


13.3 


50 


1088 


2.9 


93 


309 


8.9 


14,7,9 


UGC 128 


00 13 50.9 


+35 59 39 


24.2 


60 


57 


4509 


13.5 


133 


62 


15.5 


8,3,3 


UGC 191 


00 20 05.2 


+10 52 48 


22.7 


17.6 


39 


1139 


2.4 


97 


156 


7.7 


14,13,12 


UGC 1551 


02 03 37.5 


+24 04 32 


22.5 


20.2 


63 


2663 


3.8 


83 


111 


9.7 


5,11,12 


NGC 959 d 


02 32 23.9 


+35 29 41 


21. 9 e 


7.8 


51 


590 


1.6 


77 


64 


8.7 


6,7,7 



Note. — Units of right ascension are hours, minutes, and seconds, units of declination arc degrees, arcminutcs, and arcscconds. Col.(l): Galaxy 
name. Col. (2): Right Ascension. Col. (3): Declination. Col. (4): Central surface brightness in i?-band. Col. (5): Distance. Col. (6): Inclination. 
Col. (7): Heliocentric systemic velocity. Col. (8): Maximum radius of the DcnsePak rotation curve. Col. (9): Maximum velocity of the DcnscPak 
rotation curve. Col. (10): Position angle of major axis; see §3 for details. Col. (11): Velocity dispersion of the DcnscPak data. Col. (12): References 
for surface brightness, distance, and inclination. 
a Converted from R band assuming B ~ R — 0.9. 
b NGC 7137 = UGC 11815. 

c Converted from V band assuming B — V — 0.57. 
d NGC 959 = UGC 2002. 

e Converted from / band assuming V — I — 0.90 and B — V — 0.53. 

REFERENCES - (1) Baggctt ct al. 1998. (2) de Blok & Bosma 2002. (3) do Blok & McGaugh 1996. (4) de Blok, McGaugh, & Rubin 2001. (5) 
dc Jong 1996. (6) Heraudeau &: Simien 1996. (7) James ct al. 2004. (8) McGaugh 2005. (9) McGaugh, Rubin, & de Blok 2001. (10) Rosenberg & 
Schneider 2003. (11) Swaters & Balcells 2002. (12) Tully 1988. (13) van Zee & Haynes 2006. (14) van Zee et al. 1997. 



R0TCUR. We used R0TCUR to determine the position an- 
gle of the major axis, using published long-slit values as 
the initial guess. If the position angle could not be well- 
constrained by R0TCUR, then it was fixed to the long- 
slit value. The position angles of UGC 4325, DDO 64, 
and F583-1 were not well-constrained by R OTCUR, and 
remained fixed at the values used in IK06L The po- 
sition angle of the major axis was well-constrained by 
ROTCUR for all of the galaxies in the new sample except 
for UGC 128. In the case of UGC 128, the position an- 
gle i_was_s^tJo_t4ie_rJositiojiangle of the Hi velocity field 
of | van der Hulst et al.l (|1993f ). For reasons described in 
IK06L we did not impose a minimum error on the rotation 
curve points. However, we added in quadrature to the 
error bars on the final ROTCUR rotation curve the veloc- 
ity error from Gaussian centroiding accuracy, typically 
~1.5 km s _1 , corrected for inclination. We have also not 
corrected the rotation velocities for asymmetric drift, as 
the corrections are typ ically only ^2 km s _1 (see also 
Ide Blok fc Bosm£]l2002h . 

4. OBSERVED VELOCITY FIELDS AND DERIVED 
ROTATION CURVES 

In this section, we present the DensePak fiber posi- 
tions, observed velocity fields, and rotation curves in Fig- 
ures 1-3. A description of each galaxy is given, and the 
properties of the galaxies for which a rotation curve was 
derived are listed in Table 1. 

4.1. Extended Observations of Previously Observed 
Galaxies 

Additional DcnscP ak ob servations were made of three 
galaxies published in lK06l . In this section, we describe 
the positions of the new DensePak pointings, present 
the augmented velocity fields, and discuss the resulting 
changes in the rotation curves. 

UGC 4325— There were four new DensePak point- 
ings roughly through the galaxy center from SE to NW. 
All eight pointings are shown on the Ha image (Figure 1) 



of the galaxy. The fiber velocities were the average of the 
Ha and [Sn]A6717 lines. The position angle of the ma- 
jor axis remained fixed at the average of the position an- 
gles of previous long -slit observations (|de Blok &: Bosmal 
I2002t ISwaters et atll2003al ). With the addition of the 
new DensePak pointings, the DensePak rotation curve 
extends to ~ 60" and remains in e xcellent agreement 
with t he long-slit Ha rotation curve of lde Blok &: Bosmal 
(|200l . 

DDO 64 There were two new DensePak pointings 
along the SE end of DDO 64. All five pointings are shown 
on the Ha image (Figure 1) of the galaxy. The fiber ve- 
locities were the average of the Ha and [Sn]A6717 lines. 
The posi tion angle remain e d fixe d at the value deter- 
mined in Ide Blok &: Bosmal (|2002l ). With the addition of 
the new DensePak pointings and an Ha image, the po- 
sitions of the DensePak fibers on the galaxy have been 
updated and the center of the velocity field has moved 
13"E and 2". 5 S. The new DensePak rotation curve is 
now consistent wi t h the long-slit Ha rotation curve of 
Ide Blok fc Bosmal (|2002l ) in the inner 20". Beyond 20", 
the DensePak rotation curve shows a more linear rise 
with less scatter than its previous version. 

F 583-1— There were two new DensePak pointings, 
one to the north and one to the east, on this galaxy. All 
three pointings are sho wn on the -R-band im age (Figure 
1) of the galaxy from Ide Blok et all (|1996l ). The fiber 
velocities were the average of the Ha, [Sn]A6717, and 
[Sn]A6731 lines. The receding side of the galaxy and 
more of the minor axis are now covered with the ad- 
ditional pointings. The position ang l e rem ained fixed 
at the value listed in lMcGaugh et all (|2001l ). The new 
DensePak rotation curve shows less scatter than the 
curve derived from a single DensePak pointing, and it dis- 
plays a higher le vel of consistency with the long-slit Ha 
rotation curve of lde Blok. McGaugh. fc Rubinl (|2001f ). 



4.2. New Observations 
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Fig. 1. — Additional observations of UGC 4325, DDO 64, and F583-1. Top row: Position of DensePak array on the Ho images of UGC 
4325 and DDO 64 and the _R-band image of F583-1. Second row: Observed DensePak velocity field with new pointings. Empty fibers are 
those without detections. Third row: Updated DensePak rotation curves. Bottom row: Updated DensePak rotation curves plotted with 
long-slit Ho and Hi rotation curves. [See the electronic edition of the Journal for a color version of this figure.] 



Of the 14 galaxies observed, there were eight galaxies 
for which meaningful velocity fields could not be con- 
structed. 5 The Ha emission in these galaxies was too 
faint to be detected and/or not spread out enough across 
the fiber array. Velocity fields and rotation curves were 
derived for the remaining six galaxies in the sample, and 
each is described below. 

NGC 1131— There were three DensePak pointings 
for this galaxy. The pointings are shown on the Ha im- 
age (Figure 2) of the galaxy. Spiral arms are clearly 
visible in this galaxy. The fiber velocities were the av- 
erage of the Ha, [Nn]A6584, [Sn]A6717, and [Sn]A6731 

5 F469-2, UGC 2034, UGC 2053, UGC 11944, UGC 12048/NGC 
7292, UGC 12082, UGC 12212, and UGC 12632. 



lines. Ha emission was abundant and the majority of 
fibers detec ted emission . The inclination was fixed to 
the value in iTullv (198cf). The position angle of the ma- 
jor axis was well-constrained by RDTCUR. The rotation 
curve rises steeply out to roughly 10" then dips slightly 
before rising again. 

UGC 11820- This galaxy is both large and diffuse 
on the sky. There is a central concentration of Ha emis- 
sion, which is perhaps a bar, that runs roughly NE-SW 
and also two large, diffuse arms that extend from the end 
of the bar. There were three DensePak pointings along 
the central feature. The pointings are shown on the Ha 
image (Figure 2) of the galaxy. The fiber velocities were 
the average of the Ha, [Sn]A6717, and [Sn]A6731 lines. 
The Ha emission was sparse and only roughly half of the 
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Fig. 2.— Results for NGC 7137, UGC 11820, & UGC 128. Top row: Position of DcnsePak array on the Ha images of the galaxies. 
Middle row: Observed DensePak velocity field. Empty fibers are those without detections. Bottom row: DensePak rotation curves. The 
UGC 128 DensePak points at large r adii with a > 4 km s —1 are plotted as smaller points. The DensePak rotation curve of UGC 128 is 
plotted with the Hi rotation curve of Vcrhciien & dc Blok (1999) (stars). [See the electronic edition of the Journal for a color version of 
this figure.] 



fibers detected e mission. The inc l inatio n was fixed to 
the value listed in McGaugh ct al. (2001). The position 
angle of the major axis was well-constrained by R0TCUR. 

UGC 128 — There were three DensePak pointings 
across this galaxy. The pointings are shown on the 
Ha image (Figure 2) of the galaxy. The fiber velocities 
were the average of the Ha, [Nn]A6584, and [Sn]A6717 
lines. The Ha emission in this galaxy was sparse, but the 
emission that was present was scattered across the three 
DensePak pointings such that both the approaching and 
receding sides of the velocity field were mapped. The po- 
sition angle was fixed to the position an gle of the Hi ve- 
locity field of I van der Hulst et al.l (119931) and the inclina- 
tion was fixed to the value listed in Ide Blok fe McGauehl 
(1996). The DensePa k rotation curve is plotted w ith the 
Hi rotation curve of IVerheiien fc de Blok] (|1999f ). The 
DensePak rotation curve is largely consistent with if not 
slightly steeper than the Hi curve and does not go out 
far enough to show a clear turn-over. 

UGC 191 There were two DensePak pointings for 
this galaxy and they are shown on the Ha image (Fig- 
ure 3). The fiber velocities were the average of the Ha, 
[Sll]A6717, and [Sn]A6731 lines. There is ample Ha 



emission in the galaxy and almost every fiber had a de- 
tectio n. The inclination was fixed to the value in iTullvl 
(1988). The position angle of the major axis was well- 
constrained by R0TCUR. The rotation curve rises linearly 
and has no clear turn-over. 

UGC 1551 There were three DensePak pointings 
for this galaxy. The pointings are shown on the Ha image 
(Figure 3) of the galaxy. The fiber velocities were the av- 
erage of the Ha, [Nn]A6584, [Sn]A6717, and [Sn]A6731 
lines. Ha emission was detected in nearly a ll of th e fibers . 
The inclination was fixed to the value in ITullvl (|1988T> . 
The position angle of the major axis was well-constrained 
by R0TCUR. There is a twist in the velocity field that is 
suggestive of the presence of a bar. Noncircular motions 
are probably important inside of 15" where the rotation 
curve is mostly flat. Beyond 15", there is a linear rise in 
the rotation curve. 

NGC 959 There were three DensePak pointings 
across the length of this galaxy. The pointings are 
shown on the Ha image (Figure 3) of the galaxy. The 
fiber velocities were the average of the Ha, [Nn]A6584, 
[Sn]A6717, and [Sn]A6731 lines. Ha emission was abun- 
dant and the majority of fibers detected emission. The 
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Fig. 3. — Results for UGC 191, UGC 1551, & NGC 959. Top row: Position of DensePak array on the Ha images of the galaxies. 
Middle row: Observed DensePak velocity field. Empty fibers are those without detections. Bottom row: DensePak rotation curves. The 
open points in the UGC 191 and UGC 1551 rotation curves were excluded from the halo fits. [See the electronic edition of the Journal for 
a color version of this figure.] 



inclina tion was fixed to the value listed in iJames et al.l 
(2004). The position angle of the major axis was 
well-constrained by ROTCUR. The rotation curve is well- 
behaved with a steady rise and a turnover to Vfi at ~ 80 
km s . 

5. ZERO DISK HALO FITS 

In this section, we present the pseudoisothermal and 
NFW halo fits to the DensePak rotation curves in the 
zero disk case. By ignoring the velocity contribution from 
the baryons and attributing all rotation to dark matter, 
we are able to put an upper limit on the slope and/or con- 
centration of the halo density profile. For those galaxies 
with photometry, halo fits for three assumptions about 
the stellar mass-to-light ratio are presented in § 6. 

5.1. Halo Models 

The cuspy NFW halo and the cored pseudoisothermal 
halo are two of the most well-known competing descrip- 
tions of dark matter halos. We provide a brief description 
of each below. 

5.1.1. NFW Profile 



Numerical simulations show that the density of CDM 
halos rises steeply toward the halo center. The exact 
value of the inner slope of the CD M halo varies s light! 
depending on the simul ation (e.g. Navarro et al.l 1199 
19971; [Moore et al.lll999tlReed et al.ll2003l ; iNavarro et all 
20041 : [Diemand et al.ll2005j ). From an observational per- 
spective, there is very little to distinguish the various 
flavors of cuspy CDM halos, and we choose to fit the 
data with the NFW halo. 

The NFW mass-density distribution is described as 



Pnfw{R) — 



{R/R s ){l + R/R s ) 2 



fl 



in which pi is related to the density of the universe at the 
time of halo collapse, and R s is the characteristic radius 
of the halo. The NFW rotation curve is given by 



V(R) = V* 



'ln(l 



cx/(l 



cx) 



;[ln(l 



-c/(l 



(2) 



with x = i?/i?2oo- The rotation curve is parameterized 
by a radius i?200 an d a concentration parameter c = 
R200/RS, both of which are directly related to R s and 
pi. Here -R200 is the radius at which the density contrast 
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Fig. 4. — Halo fits to the DensePak rotation curves. The red solid line is the best-fit isothermal halo, the green short-dashed line is the 
best-fit unconstrained NFW halo, and the black long-dashed line is the best-fit NFW CO nstr halo. NFW con str fits were only made to UGC 
4325, DDO 64, and F583-1. [See the electronic edition of the Journal for a color version of this figure.] 



exceeds 200, rou ghly the virial radius: V2no is the circular 
velocity at R 20 o (|Navarro et al.lll996lll997h . Because the 
NFW profile has a shallower slope than other cuspy halo 
models, it provides a lower limit on the slope of cuspy 
density profiles and, as such, gives the cuspy halo the 
best possible chance to fit the data. 

5.1.2. Pseudoisothermal Halo 

The pseudoisothermal halo describes a dark matter 
halo that has a core of roughly constant density. By con- 
struction, it produces flat rotation curves at large radii. 
The density profile of the pseudoisothermal halo is 

p lso (R) = Po [l + (R/Rc) 2 r\ (3) 
with po being the central density of the halo and Rc 
representing the core radius of the halo. The rotation 
curve corresponding to this density profile is 



V(R) 



UirGpoR 2 



c 



1 



(4) 



Rc , ( R 

— — arctan — — 

R \R C/ 

The pseudoisothermal halo is empirically motivated and 
predates halo profiles stemming from numerical simula- 
tions. 



5.2. Halo Fits to Previously Observed Galaxies 

We find the best-fit zero disk case isothermal and NFW 
halos to the new DensePak rotation curves of UGC 4325, 
DDO 64, and F583-1. When available, the DensePak 
rotation curves have been supplemented with previous 
smoothed long-slit Ha and Hi rotation curves. We use 
the entire long-slit rotation curve, and include only those 
Hi points that extend beyond the radial range of both the 
DensePak and long-slit data. Uncertainties from possible 
resolution effects are avoided by using only the outer Hi 
points. 

For these three galaxies, we also fit an NFW halo called 
NFW mr istTO i n ed (hereafter NFW COIlstr ), as described in 
detail in lKOtj . This halo fit is motivated by NFW fits that 
have parameters which are unrealistic or inconsistent 
with ACDM. Briefly, the constrained halo was required 
to match the velocities at the outer radii of each galaxy 
while constraining the concentration to agree with cos- 
mology. T he concentrati o ns we re calculated using Equa- 
tion (7) of Ide Blok et~ai] (j200j i wlrich gives the co ncen- 
tration as a function of of V200 ilNayarro et al.lll997), and 
then adjusted to the cosmology of iTegmark et alJ (2004) 
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TABLE 2 
Isothermal Halo Parameters 



Galaxy 




ZERO DlbK 








MINIMAL DlbK 






R c 


Pa 


Xr 


T* 


R c 


Pa 


xt 


T, 




4.1±0.3 


gg-j-3 


3 2 





4.6±1.5 


77±5 


3 1 


0.57 


r juo v ^ 


1.5±0.1 


119±6 


0.71 





1.4±0.2 


114±20 


0.68 


0.44 


F563-1 


21-J-0.1 


67±2 


0.43 





2.0±0.2 


61il0 


0.48 


0.68 




3 3±0.5 


43±3 


3 2 





4.1±3.8 


37±6 


3.3 


0.62 


F568-3 


3.8±0.2 


27±1 


1.2 


0.0 


4.0±0.5 


22±3 


1.5 


0.66 


UGC 5750 


5.7±0.4 


7.1±0.3 


0.83 


0.0 


6.5±1.1 


5.2±0.6 


0.84 


0.62 


NGC 4395 


0.7±0.1 


258±9 


2.9 


0.0 


0.57±0.05 


318±42 


2.7 


0.70 


F583-4 


1.3±0.1 


67±2 


0.67 


0.0 


1.2±0.2 


66±16 


0.62 


0.53 


F583-1 


2.5±0.1 


30±2 


0.50 


0.0 


2.4±0.2 


30±3 


0.59 


0.62 


NGC 7137 


0.6±0.1 


274±17 


3.1 


0.0 










UGC 11820 


l.liO.l 


274±21 


2.9 


0.0 










UGC 128 


2.3±0.1 


66±0.6 


8.5 


0.0 










UGC 191 


1.7±0.1 


138±10 


5.8 


0.0 










UGC 1551 


1.3±0.2 


57±5 


5.5 


0.0 










NGC 959 


0.4±0.1 


1117±29 


1.2 


0.0 











POPSYNTH MAXIMUM DISK 

Galaxy R c pp \ r X* Rc P" *r T * 



UGC 4325 a 6.2±3.3 66±5 3.2 1.14 2.5 45 11 4.5 

F563-V2 1.4±0.2 102±20 0.75 0.88 4.7±5.9 7.1±8.8 1.4 4.0 

F563-1 6 2.1±0.2 54±10 0.50 1.36 4.2±1.2 13±6 0.83 6.5 

12.5±9.0 2.7±1.6 1.3 10.0 

DDO 64 a 5.9±11.1 32±5 3.5 1.24 1.5 24 2.2 5.0 

F568-3 4.7±0.8 16±3 1.8 1.32 6.1±1.8 10±3 2.5 2.3 

UGC 5750 7.5±1.6 4.1±0.5 0.95 1.24 9.8±3.2 2.6±0.5 1.3 2.2 

NGC 4395° 0.50±0.04 355±50 2.8 1.40 19 0.32 4.1 9.0 

F583-4 a 1.2±0.2 63±16 0.62 1.06 8 1.4 1.2 10.0 

F583-1 6 2.5±0.2 27±3 0.60 1.24 3.6±0.5 14±2 0.83 5.0 

6.9±1.7 5.2±0.9 2.0 10.0 



Note. — R a is in kpc; p is in 10~ 3 M Q pc -3 . Photometry is unavailable for NGC 7137, UGC 11820, UGC 128, UGC 191, UGC 1551, & NGC 
959; halo fits beyond zero disk arc not presented. 

a The baryons can account for most of the velocity in the maximum disk fit. See text for details of fit. 
As discussed in the text, F563-1 and F583-1 have two possible values of T^Max- 



by su btracting 0.011 dex (|McGaugh, Barker, fc de Blokl 
12001 . 

The halo fits are plotted over the data in Figure 4, and 
the halo parameters are listed in Tables 2 and 3. For 
comparison, th e num bers mentioned in the text below 
are the fits from lK06l . unless specifically noted otherwise. 

UGC 4325- The isothermal fit to UGC 4325 im- 
proves slightly with the addition of the four new 
DensePak pointings {R c = 3.3 ± 0.2; p Q = 91 ± 4; X 2 
= 3.8). An unconstrained NFW halo still could not be 
fitted to the data, and the poor quality of the NFW co „ str 
fit remains virtually unimproved (c = 6.9; V200 = 249; 
X 2 = 40). This galaxy remains best-described by the 
isothermal halo. 

DDO 64 — The two new DensePak pointings on DDO 
64 help to improve the quality of both the isothermal and 
NFW con;str fits. The values of the isothermal halo param- 
eters remain within the errors of the original values, but 
the value of Xr decreases (R c — 4.4 ± 0.9; po = 38 ± 3; 
X 2 = 5.5). While no unconstrained NFW halo fit could 
be made, the quality of the NFW CO nstr halo significantly 
improves (c = 9.2; V200 = 62; x 2 = 20). This does not 
mean, however, that the NFW halo is a good fit to the 
data. The NFW constr halo continues to overshoot the 
data at radii interior to where it was forced to match the 
data. While it passes through the large errorbars on the 
inner rotation curve points, significant noncircular mo- 
tions would need to be important all the way out to ^1 
kpc in order to boost the observed velocities up to the 



expected NFW velocities. 

F 583-1- The isothermal and NFW fits are signifi- 
cantly better constrained by the addition of the two new 
DensePak pointings. F583-1 remains best-described by 
the isothermal halo. The values of the halo parameters 
remain within the errors of the previous values, but the 
value of Xr falls from 5.4 to 0.5 (R c = 2.7 ± 0.1; p = 35 
± 2). The new values of the isothermal halo parameters 
ar e almost indist i nguish able from the values determined 
bv lde Blok et~aTl (pOOl : R c = 2.44 ± 0.06, p a = 33.0 ± 
1.1. The values of the parameters of the unconstrained 
NFW fit are also mostly unchanged, and Xr decreases to 
1.7 from 8.7 (c = 4.7 ± 0.7; V 20 a = 133 ± 21). The value 
of the concentration, c = 4.5, is still too low to be con- 
sistent with ACDM. The value of x 2 a ls° drops for the 
NFW constr fit (x 2 = 11 -> 6.5), but the NFW velocities 
continue to over-predict the observed velocities interior 
to where they were forced to match the data. 

5.3. Halo Fits to New Observations 

We find the best-fit zero disk case isothermal and 
NFW halos to the DensePak rotation curves of NGC 
7137, UGC 11820, UGC 128, UGC 191, UGC 1551, and 
NGC 959. In the case of UGC 128, we combine the 
DensePak rotation curve w ith the Hi rotation curve of 
IVerheijen fc de Blokl (fl999l. using only those Hi points 
beyond the radial range of the DensePak data. We do 
not make an NFW cons t r fit to these six galaxies, as the 
radial range of the data does not extend beyond the ris- 
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TABLE 3 
NFW Halo Parameters 







ZbKU DISK 








MINIMAL DISK 




Galaxy 


c 


V200 


— — 


Y* 


c 


^200 


Xr 


Y„ 


HHP 4395 


6 9 


249 


39 





1 


1002 


12 


0.57 


IT UUO V ^ 


7.7±2.0 




0.40 





8.4il.5 


105±17 


0.46 


0.44 


F563-1 


7.8±1.3 


106il0 


0.88 





7.6il.3 


100±9 


0.89 


0.68 


DDO 64 


9.2 


62 


12 


0.0 


1.0 


376 


6.5 


0.62 


F568-3 


8.2 


110 


12 


0.0 


1.0 


465 


3.9 


0.66 


UGC 5750 


0.5±0.1 


320±43 


1.7 


0.0 


1.0 


167 


1.7 


0.62 


NGC 4395 


10.1±0.6 


77±4 


2.1 


0.0 


11.5±1.0 


63±4 


2.1 


0.70 


F583-4 


5.5±2.2 


92±32 


0.41 


0.0 


5.7±1.4 


83±18 


0.41 


0.53 


F583-l a 


4.5±0.8 


120±20 


1.7 


0.0 


5.1±1.2 


102±21 


1.8 


0.62 


NGC 7137 


15±3 


56±10 


3. 1 


0.0 










UGC 11820 


















UGC 128 


8.9±0.2 


111±0.7 


9.3 


0.0 










UGC 191 


















UGC 1551 


















NGC 959 


23±4 


76±15 


1.7 


0.0 
















POPSYNTH 








MAXIMUM DISK 




Galaxy 


c 


V200 


xi 


T» 


c 


V200 


Xr 


T* 


UGC 4325 


1.0 


897 


12 


1.14 


1.0 


283 


12 


4.5 


F563-V2 


7.8±1.7 


104±20 


0.52 


0.88 


1.0 


203 


1.1 


4.0 


F563-1 6 


7.0±1.3 


102±11 


0.90 


1.36 


1.0 


283 


1.0 


6.5 












1.0 


210 


1.5 


10.0 


DDO 64 


1.0 


332 


6.4 


1.24 


1.0 


146 


2.6 


5.0 


F568-3 


1.0 


404 


4.3 


1.32 


1.0 


317 


5.2 


2.3 


UGC 5750 


1.0 


141 


1.9 


1.24 


1.0 


101 


2.4 


2.2 


NGC 4395 


12.5±1.1 


61±4 


2.2 


1.40 








9.0 


F583-4 


5.8±1.5 


86±19 


0.40 


1.06 








10.0 


F583-1 6 


4.9±1.2 


110±26 


1.9 


1.24 


1.0 


298 


2.1 


5.0 












1.0 


173 


3.3 


10.0 


Note. — V 2 oo 


is in km s — 1 . 


Italicized halo parameters are 


forced fits. See 


text for details. 


Photometry is unavailable for NGC 7137, UGC 11820, 



UGC 128, UGC 191, UGC 1551, & NGC 959; halo fits beyond zero disk are not presented. 

a The parameters of the updated N FWcanstr fit for the zero disk case of F583-1 are c — 8.7, V200 — 83, Xr — 6.5. 
b As discussed in the text, F563-1 and F583-1 have two possible values of Y* May- 



ing part of the rotation curve. In Figure 4 we plot the 
halo fits over the data, and list the halo parameters in 
Tables 2 and 3. 

NGC 1131- The isothermal halo is a slightly better 
fit to NGC 7137 than the NFW halo, but the difference 
is not particularly significant. There are "bumps and 
wiggles" in the rotation curve that simple, smooth halo 
models cannot fit. The value of the concentration, c = 
15, for the NFW fit is on the high side of values expected 
in a ACDM cosmology. 

UGC 11820- UGC 11820 is fit relatively well by the 
isothermal halo; no NFW halo could be fit to the rotation 
curve. 

UGC 128- UGC 128 is more consistent with the 
isothermal halo than the NFW halo, though the value 
of Xr is high for both fits and both halo models overes- 
timate the rotation velocities at small radii. The best- 
fitting concentration, c = 8.9, is reasonable for a galaxy of 
this size in a ACDM cosmology. The best-fitting isother- 
mal halo parameters for the combined DensePak+Hi 
rotation curve are R 9 = 2.3±0.1 and po — 66±0.6. 
de Blok & McGaugh (1996) also fit an isothermal halo 
to the Ivan der Hulst et al.l ((1993) Hi rotation curve and 
find R c = 4.0 and po = 21.7. The differences in the halo 
parameters reflect the inclusion of the slightly steeper 
DensePak rotation curve. 

UGC 191 We exclude the innermost point of the 
UGC 191 rotation curve due to its unrealistically small 
formal error bar and the steep jump in velocity between 



it and the next rotation curve point. Excluding this point 
from the halo fits does not significantly alter the values of 
the halo parameters, but does improve the values of the 
reduced \ 2 ■ No NFW fit could be made to the DensePak 
rotation curve. The data were fit with an isothermal 
halo, albeit with a large \r- The UGC 191 DensePak 
rotation curve is linearly rising and shows no turn-over 
over in the radial range covered by the DensePak data. 
To be useful for distinguishing between halo models, data 
points at larger radii are needed to better constrain the 
halo fits. 

UGC 1551 We exclude the first five points of the 
UGC 1551 rotation curve because the twist in the ve- 
locity field suggests that a bar may be present in the 
galaxy. An NFW fit could not be made to the data; an 
isothermal fit was made, but was not well-constrained. 
Like UGC 191, data at larger radii are necessary for ob- 
taining useful constraints on the halo models. 

NGC 959- NGC 959 is well described by the isother- 
mal halo and is a slightly better fit than the NFW halo. 
The value of the concentration, c = 23, of the best fit 
NFW halo is on the high side of expected values for 
galaxies in a ACDM cosmology. 
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5.4. Summary of Zero Disk Halo Fits 

In § 5.2 and § 5.3 we have pr esent ed updated zero disk 
halo fits to three galaxies from IK06I and six new galaxies 
observed with DensePak. Overall, we find the isothermal 
halo to be a better description of seven of these galaxies 
than the NFW halo. Two of the new galaxies do not 
have data at large enough radii to put useful constraints 
on the halo models. When NFW fits could be made, 
the concentrations were often beyond the range of val- 
ues expected for ACDM. We find that the parameters of 
the halo fits to the three previously observed galaxies do 
not significantly change, but are better constrained by 
the additional DensePak coverage. The DensePak rota- 
tion curves of the six new galaxies could all be fit by the 
isothermal halo, and an NFW fit could be made to only 
three. Only one of the three NFW fits had a concentra- 
tion consistent with the range expected for galaxies in a 
ACDM cosmology; the other two are on the high end of 
expected values. We also find that the quality of both 
the isothermal and NFW halo fits is greatly improved 
when the radial range of the data extends into the flat 
part of the rotation curve. Of the six new galaxies, only 
NGC 959 has a DensePak rotation curve which clearly 
turns over and flattens, and it is this galaxy that has the 
isothermal and NFW fits with the lowest Xr- 

6. MASS MODELS 

Low surface brightness galaxies are dark matter- 
dominated down to small radii. Because of this, the ve- 
locity contribution from the baryons is often disregarded 
when fitting dark matter halo models to the galaxy rota- 
tion curves. This type of fit which ignores the contribu- 
tion of the stars and gas, such as in § 5, is usually called 
the minimum-disk case. While the dark matter is the 
dominant mass component of LSB galaxies at all radii, 
baryons are still important. To accurately determine the 
distribution of the dark matter, it is necessary to prop- 
erly account for the stars and gas in the galaxies and 
their contribution to the observed rotation. The veloc- 
ity component coming from the stars is determined from 
the surface photometry scaled by the stellar mass-to-light 
ratio, T*. This is a straightforward computation of the 
gravitational potential of the observed stars. However, 
the true value of Y* is difficult to determine, and as dis- 
cussed below, can be assigned a value following a number 
of techniques. Similarly, the velocity component from the 
gas is determined using Hi surface density profiles. The 
conversion from observed 21 cm luminosity to atomic gas 
mass is well understood from the physics of the spin-flip 
transition. A scaling factor, the inverse of the hydro- 
gen mass fraction, is usually included to account for the 
helium and metals also pres ent in the galaxies. 

In § 5 of this paper and in lK06l . we have so far ignored 
the baryons in the halo fits to LSB galaxies observed with 
DensePak. For those galaxies with i?-band photometry 
and Hi surface density profiles we now include the con- 
tribution from the baryons and present the dark matter 
halo fits for three determinations of T*. 

6.1. Dynamical Components 

There are three separate components in galaxy mass 
models that contribute to the observed velocities: the 
stars, the gas, and the dark matter. The stars, gas, and 



TABLE 4 

Galaxy Parameters 



Galaxy 


h (kpc) 


(B-R) 


r» (Pop) 


References 


UGC 4325 


1.6 


0.85 


1.14 


1,3 


F563-V2 


2.1 


0.51 a 


0.88 


2,4 


F563-1 


2.8 


0.96 


1.36 


2,5 


DDO 64 


1.2 


0.9 6 


1.24 


1 


F568-3 


4.0 


0.94 


1.32 


2.5 


UGC 5750 


5.(5 


0.9 b 


1.24 


2.2 


NGC 4395 


2.3 


C 


1.40 


1 


F583-4 


2.7 


0.8 


1.06 


2.6 


F583-1 


1.6 


0.9 


1.24 


2.6 



Note. — Col. (2): Scale length (kpc). Col. (4): T. for the popsynth 
case as determined from the colors. Col. (5): References for h and (B- 
R), respectively: (1) dc Blok & Bosma (2002) (2) do Blok, McGaugh, 
& Rubin (2001) (3) van den Bosch & Swatcrs (2001) (4) de Blok & 
McGaugh (1997) (5) dc Blok, van dor Hulst, & Bothun (1995) (6) dc 
Blok, McGaugh, & van der Hulst (1996). 
a (B-V) color. 

b Multicolor photometry unavailable; assuming (B-R)— 0.9 for dwarf 
galaxies. 

c Multicolor photomety unavailable; assuming T* — 1.4. 

dark matter are added together in quadrature to obtain 
the total velocity, V t 2 otal = T*vl + V g 2 as + V* M . In this 
section, we describe each of these dynamical components. 

6.1.1. Stars 

Galaxy photometry in combination with a stellar mass- 
to-light ratio, T„, is used to determine the stellar con- 
tribution to the observed galaxy rotation. It is difficult 
to determine the true value of T» of a galaxy because it 
depends on many factors ranging from the initial mass 
function to the extinction. To address this uncertainty, 
we consider four different scenarios for the value of T*. 
We assume that T* is constant with radius. 

Zero disk— In this limiting case we ignore the con- 
tribution of the stars (T*= 0) and gas and attribute all 
rotation to dark matter. This case is typically called 
'minimum disk' in the literature; we have chosen to refer 
to this case as zero disk because it aptly describes the 
omission of the baryons. The zero disk case puts an up- 
per limit on the slope and/or concentration of the halo 
density profile. The res ults fo r the zero disk case were 
presented in § 5 and also lK06l (where they were referred 
to as minimum-disk) and are reproduce d in T ables 2 and 
3 for completeness. Two galaxies from lK06l UGC 1281 
and UGC 477, are excluded from mass models beyond 
zero disk because of their high inclinations and possible 
associated line-of-sight integration effects. Photometry 
is unavailable for the six new galaxies presented in § 5; 
mass models beyond zero disk will not be presented. 

Minimal disk— Though commonly employed, the 
zero disk case is unphysical. We next consider a more 
realistic minimal contribution of the stars and the gas. 
The Y» is determined from population synthesis models 
(see below) and then scaled by 0.5 (see lMcGaugh|[2005l) 
to simulate a lightweight IMF. 

Popsynth Here T» i s determined usin g the popula- 
tion synthesis models of iBell et al.l (|2003h . Specifically, 
we use the relation between T* and color as defined in 
their Table 7. For a scaled Salpeter IMF, (B-R) and 
(B—V) colors are related to T*. This is our best esti- 
mate of the baryonic mass from the perspective of stellar 
populations. The popsynth model also includes the gas 
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Fig. 5. — Isothermal and NFW halo fits for UGC 4325. The green dotted line is the rotation curve of the gas disk, the gold short-dashed 
line is the rotation curve of the stellar disk, the blue long-dashed line is the rotation curve of the dark matter halo, and the red solid line 
is the total model curve. The velocity contribution from both the stars and gas is ignored in the zero disk case. The minimal disk case 
considers the gas contribution and assumes a lightweight IMF for the stars. The popsynth case includes the gas contribution and uses 
population synthesis models to determine the stellar contribution. The maximum disk case includes the gas contribution and the stellar 
mass-to-light ratio is scaled up as far as the data will allow. [See the electronic edition of the Journal for a color version of this figure.] 



contribution. Colors and corresponding T* are listed in 
Table 4. 

Maximum disk— We also consider the case where 
the T* is scaled up as far as the data will allow. In high 
surface brightness galaxies, this approach usually fits the 
inner rotation curve well. In LSB galaxies, the shape of 
the stellar rotation curve is often not well-matched to the 
observed rotation curve. We have chosen a that al- 
lows the stellar rotation curve to match the inner rotation 
curve points as well as possible, sometimes overshooting 
the innermost point in order to hit the next few. These 
Y* are usually higher than those determined from the 
popsynth models, and in some cases (eg. F583-4) are 
substantially higher. 

To model the stellar disk we have u s ed the J?-band pho - 
tometry presented in Ide Blok et alj (|1995t ). IStU (l999), 



iSwaterd £1999), and iSwaters et all (|2002f ). The GIPSY 
task ROTMOD was used to determine the rotation of 
the disk assuming a ve rtical sech 2 distribution wit h 
a scale height ZQ—h/6 (jvan der Kruit fc Searld fl98lh . 
The stellar rotation curve computed from the photom- 
etry was resampled at the same radii as the combined 
DensePak+long-slit+Hi rotation curves. 

Each estimator of Y* has its advantages and disadvan- 
tages. The zero disk T* is often adopted and is useful 
because it provides an upper limit on the dark matter. 
Unfortunately, it is an unphysical assumption. Popu- 
lation synthesis models best represent what we know 
about stars, however, sometimes the data will allow a 
larger Y*. And although T^max often seems too large 
with respect to Y*pop, disk featu res like bars and spi- 
ral arms require a large disk mass (|McGaugh fc de B lok 
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Il998t lFuchsll2003fl and support high T *MAX- 
6.1.2. Gas 

The gas present in galaxies also contributes to the ob- 
served galaxy rotation. Hi is the dominant gas com- 
ponent, but to include helium and metals, the Hi data 
were scaled by a factor of 1.4. Substantial amounts of 
molecular gas are not obviously present in LSB g alaxies 
(|de Blok fc van der Hulsil 119981 : ISchombert et al.lll990D . 
The total m ass of H2 is almost certainly much less than 
that of Hi (jMihos. Spaans. fc McGaughl I1999D . More- 
over, in brighter ga laxies H2 is known to trace the stars 
(|Regan et al.ll200lh . so at most represents a slight tweak 
to T». The contribution of the gas to the observed ve- 
locity is considered in the minimal, popsynth, and max- 
imum disk cases. 

The Hi s ur face den sity profiles pres ent ed i n 
Ide Blok et all (11996D . Ivan der Hulst et all (|1993l ). 
and iSwaterj (| 19991) were used to model the gas disk. 



The GIPSY task ROTMOD was used to determine the 
rotation of the disk assuming a thin disk. The gas 
rotation curve was resampled at the same radii as the 
combined DensePak+long-slit+Hi rotation curves. 

6.1.3. Dark Matter Halo 

In LSB galaxies there is usually a considerable amount 
of observed rotation unaccounted for after subtracting 
off the velocity of the stellar and gas disks, even in the 
maximum disk case. The remaining rotation is usually 
attributed to dark matter. While there have been a num- 
ber of dark matter models proposed in the literature, we 
fit two of the most prominent competing profiles: the 
pseudo-isothermal halo and the NFW profile. Both halo 
profiles are described in § 5.1. 

7. MASS MODEL RESULTS FOR INDIVIDUAL GALAXIES 

In this section, we present the isothermal and NFW 
halo fits for the four scenarios of T*. The halo parame- 
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Fig. 7. — Isothermal and NFW halo fits for F563-1. Line types are described in Figure 5. [See the electronic edition of the Journal for 
a color version of this figure.] 



ters are listed in Tables 2 and 3. Figures 5-13 plot the 
halo fits over the data. 

As T* increases and the baryons become responsible 
for more of the observed velocity, less room is available 
for dark matter and the NFW halo becomes increasingly 
difficult to fit to the data. The best-fit concentrations 
drop to very small, and sometimes negative, values. In 
these cases, we forced an NFW fit with c = 1.0. The halo 
parameters of forced fits are italicised in the tables. 

UGC 4325- UGC 4325 is clearly best described by 
an isothermal halo; an NFW halo could not be fit to the 
data for any Y* scenario. With T^max — 4.5, the stel- 
lar rotation curve is able to describe the data well out 
to a radius of ~ 1.3 kpc, as compared to < 0.4 kpc for 
T *pop — 1-14. Because the baryons are able to explain 
so much of the data in the maximum disk case, there is 
little room left for a dark matter isothermal halo. When 
this happens, the cored halo becomes very nearly hollow 
with po decreasing to a very small value and R c increas- 



ing to a large number. For the maximum disk case of 
UGC 4325, we forced an isothermal halo fit by fixing Vh 
to approximately the observed velocities of the outer ro- 
tation curve and then varying R c such that Vtot follows 
the data. In this fit, the resulting R c is a lower limit and 
po is an upper limit. 

F563- V2— As discussed in lK06l for the zero disk case, 
F563-V2 has too few data points to really distinguish be- 
tween halo types. This remains true for the other T* sce- 
narios as well. The value of T* can be turned up to 4.0 
in the maximum disk case, leaving essentially no room 
for an NFW halo; a fit was forced with c = 1.0. NFW 
halos could be fit to the zero, minimal, and popsynth 
cases, however, and the values of the best-fitting con- 
centrations, c — 7.7, 8.4, and 7.8, respectively, are com- 
parable to values expected from simulations. Because 
of the C-V200 degeneracy that allows halos of different 
c, Vmn ) to look the same over a finite range of radius 
de Blok et all l2001f) . these three NFW fits are essen- 
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Fig. 8. — Isothermal and NFW halo fits for DDO 64. Line types are 
a color version of this figure.] 



tially indistinguishable. This is an example of where the 
zero disk assumption is reasonable for LSB galaxies. 

F 563-1 - The isothermal halo fits the F563-1 data 
better than the NFW halo, though the values of the best- 
fitting concentrations are consistent with values expected 
in ACDM. The baryons can account for the majority of 
the observed rotation in the maximum disk case and a 
forced NFW fit was made with the concentration fixed 
at 1.0. Stellar population models set T*pop = 1-36 for 
F563-1. With this value of the T„, the stellar rotation 
curve just grazes the innermost observed rotation curve 
point. The T»max can be substantially turned up such 
that the stellar rotation curve matches the data more 
closely. If T*max — 6.5, V* overshoots the innermost 
rotation curve point, but matches the cluster of points 
at ~ 1.5 kpc and even crosses the lower errorbars on 
points between 5 kpc and 10 kpc. The shape of the ob- 
served rotation curve is well- matched when T*max = 
10.0; V* goes through the upper errorbars on points < 2 
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described in Figure 5. [See the electronic edition of the Journal for 



kpc and agrees well with the data between 5 kpc and 10 
kpc. Fitting the inner observed rotation curve points or 
the overall rotation curve shape are both equally plausi- 
ble a pproaches to defining Y*max ()Palunas fc Williams] 
2000), so we use both values in our halo fits. 

F563-1 is a good example of how the choice of T* af- 
fects the mass models, and the degeneracy between the 
luminous and dark components. While T*pop is our 
best estimate for the stars, the data will clearly allow 
a higher T*. To constrain the degeneracy between the 
stars and dark matter in cases like this, we need more 
information. The stellar velocity dispersion perpendicu- 
lar to the disk for instance, would h elp to put limits on 
the disk mass (|Verheiien et al.|[2007f ). 

DDO 64 — As Hi surface density profiles were unavail- 
able for this galaxy, the baryons in the mass models are 
represented by the stars only. Multicolor photometry 
was also unavailable, so a ( B-R) = 0.9 color was as- 
sumed for this dwarf galaxy dde Blok et al.|[l995[) . DDO 
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Fig. 9. — Isothermal and NFW halo fits for F568-3. Line types are 
a color version of this figure.] 

64 is better described by the isothermal halo than the 
NFW halo; forced c — 1 NFW fits are made for each 
value of T*. The stellar rotation curve falls below the 
observed rotation curve with the exception of two low 
points in the popsynth case, T*pop = 1.24. In the max- 
imum disk case with T*max =5.0, however, V* is able 
to follow the data very well out to just past 1 kpc, and 
is even consistent with a data point at ~ 1.5 kpc. With 
TtMAX = 5.0, the stars are able to account for most of 
the observed rotation, and the displayed isothermal halo 
fit is the upper limit on po arL d the lower limit on R c . 

F568-3 — F568-3 is fit well by isothermal halos; only 
forced NFW fits could be made to the data. The shape 
of the stellar rotation curve is not well-matched to the 
observed rotation curve and is only able to describe the 
inner 2 kpc of data, even in the maximum disk case. 

UGC 5150 Because multicolor photometry was un- 
availabl e, a (B-R) = 0-9 color was assumed for this 
galaxy (|de Blok et al.l ll995). Excellent isothermal fits 
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described in Figure 5. [See the electronic edition of the Journal for 



were made for UGC 5750; the NFW halo provides very 
poor fits. Like F568-3, there is not a substantial differ- 
ence between T*pop and T*max- 

NGC 4395 Multicolor photomet ry was unavailable 
for t his galaxy and T*pop = 1.4 (de Blo k fc Bosnia! 
2002) was assumed for the popsynth model. The NFW 
halo is a slightly better fit to NGC 4395 than the isother- 
mal halo, and the best-fitting concentrations are consis- 
tent with expectations from ACDM. There is a substan- 
tial difference between T*pop and T*max- For Y»pop 
= 1.4, V* is well below the observed rotation curve; how- 
ever, for T*max — 9.0, K is able to trace the data out 
to 8 kpc, the entire length of the rotation curve. Because 
the baryons can explain the observed rotation so well in 
the maximum disk case, the data want an isothermal 
halo with an almost hollow core. We force a fit with an 
upper limit on po and a lower limit on R c . The baryons 
do such a good job of explaining the observed rotation in 
the maximum disk case that there is not even room for 
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Fig. 10. — Isothermal and NFW halo fits for UGC 5750. Line types are described in Figure 5. [See the electronic edition of the Journal 
for a color version of this figure.] 



an NFW halo with c = 1.0. In the NFW maximum disk 
plot in Figure 11, the magenta line represents the total 
baryonic rotation curve, V^ isk — V? + Vg as . It should 
be noted that the nucle us of NGC 4395 is the least lumi- 
nous known Seyfert 1 (|Filippenko &: Sargent! IT989I ) and 
may influence the inner rotation curve that we derive. 

F583-4 — The isothermal halo is a good description of 
the data. The NFW halo fits are comparable, although 
the values of the concentrations are on the low side of 
expected values from simulations. The stellar rotation 
curve is far below the data at all radii for the popsynth 
T*pop = 1-06. The entire observed rotation curve can be 
well described by the stellar rotation curve when Y*max 
is turned up to 10.0. In this case, there is very little 
room left for a dark matter halo, and only a limiting 
isothermal halo fit is made. Like NGC 4395, there is 
not enough velocity left for an NFW halo with c = 1.0, 
and the magenta line in the maximum disk NFW plot in 
Figure 12 represents the total baryonic rotation curve. 



F 583-1- F583-1 is better fit by the isothermal halo 
than the NFW halo. Additionally, the best-fitting NFW 
concentrations are on the low side of expected values 
from ACDM. Like F563-1, we consider two values of 
^*max- With T* max = 5.0, V* goes through the data 
within ~ 1.5 kpc. V* goes through the upper errorbars 
on the inner rotation curve, through the data at ~ 2 kpc, 
and then through the lower errorbars on the data out to 
~ 3.5 kpc when Y*mAX = 10.0. 

8. DISCUSSION 

While one would like to know the true T* for each 
galaxy, our data do not indicate that a particular esti- 
mator of T* is any better than another. In Figure 14 
we show that the parameters of both the isothermal and 
NFW halo fits do not change much as T* changes. This 
confirms that the details about what is assumed for the 
stars in LSB galaxies do not really matter. Unfortu- 
nately, it also means that without additional informa- 
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Fig. 11. — Isothermal and NFW halo fits for NGC 4395. Line types are described in Figure 5. The solid magenta line in the maximum 
disk NFW plot is the total baryonic rotation curve. See text for details. [See the electronic edition of the Journal for a color version of 
this figure.] 



tion, we cannot constrain T* in galaxies like F563-1 and 
F583-1 where a wide range of T* are applicable. 

Though the exact assumption about the stars may be 
unimportant, the stars cannot be entirely ignored. In 
reality, stars do not have zero mass. In fact, as the ve- 
locity contribution from the stars becomes more impor- 
tant (ie. as T* goes up), there is less room for dark 
matter at the centers of the galaxies. The left panel of 
Figure 14 illustrates that the isothermal core radius is 
often larger than ~ 1 kpc. There is also a trend towards 
larger R c with increasing Y». This is important to rec- 
ognize, as it shows that the cusp-core problem is not 
restricted to the innermost radii only, particularly when 
stars are allowed to have mass. Also plotted in this fig- 
ure is the Pq-R c scaling relation of lKormendv fc Freeman! 
(2004) (their equation 20 ) and the Pq-Rq relation de- 
termined by ISpano et~aT1 (120071 ) (their figure 2). The 
iKormendy fc Freeman! (|2004h relation follows the data 



more closely than the ISpano et a l. (2007) relation, how- 
eve r the maximum disk case poi nts tend to drift below 
the IKormendy fc Freeman! (|2004l ) relation. 

The c- V200 plot in the right panel of Figure 14 simi- 
larly shows the difficulty of fitting the centrally concen- 
trated NFW halo to the data as the velocity contribution 
from the baryons becomes larger: as T, increases, fewer 
galaxies can be fitted with an NFW halo. In particular, 
the velocity contribution from the baryons is significant 
enough in the maximum disk case that maximum disk 
and the NFW halo are mutually exclusive. 

For comparison, we have also include d the predicted 
c- V2QQ lines for the high er 8 cosmology of iTegmark et al.l 
(120041 hereafter T04 ) and the low er 8 WMAP3 cosmology 
(jSpergel et aTll2007l ). The width of the bands is ±la 
assuming a scatter of er c = 0.14 (|Bullock et alJ feOOl). 
The data appear to follo w a steeper c- V200 rela tion than 
both predictions (see also lMcGaugh et al.ll2007t ). and the 
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Fig. 12. — Isothermal and NFW halo fits for F583-4. Line types are described in Figure 5. The solid magenta line in the maximum disk 
NFW plot is the total baryonic rotation curve. See text for details. [See the electronic edition of the Journal for a color version of this 
figure.] 

concentrations of the DensePak galaxies with NFW fits 
primarily cluster between the WMAP3 lines. The me- 
dian concentration of all the zero disk DensePak galax- 
ies is c=4.5. At the corresponding V200, the median 
c expected from T04 is ~8.2, and the median c ex- 
pected from WMAP3 is ~6.1. An alternative method 
of measuring the halo central density is the A y/ 2 ap - 
proach proposed bv lAlam. Bullock. &; Weinberg) (2002). 
The median Ay/ 2 of all the zero disk DensePak galax- 
ies is ^2.9xl0 5 . At the corresponding V max , the median 
Ay/2 expected from T04 is ~4.5xl0 5 and is ~1.8xl0 5 
from WMAP3. The low concentrations/central densities 
observed in the data are more consistent with a power 
spectrum havi ng a lower amplitude on small (galaxy) 
scales (see also lMcGaugh et aLll2003|) . 

8.1. Cusp Mass Excess 



In § 5.2 we defined a constrained NFW halo, 
N FW cons tr, which was constructed to have a 
cosmologically-consistent concentration by forcing 
the halo to match the velocities at the outer radii of 
each galaxy. We can now find N FW cons t r halos for 
the dark matter rotation curves for each determination 
of Y*. The rotation curves of the NFW cons t r halos 
over-predict the observed velocities interior to where 
they are forced to agree. We can evaluate the difference 
between this expected CDM rotation curve and the 
observed dark matter rotation curve in terms of velocity 
difference, or alternatively, as a cusp mass excess. We 
can ask what the difference is between the expected 
cuspy NFW halo mass and the dark matter mass that 
is allowed by the data. In the left panel of Figure 15, 
we show the observed dark matter rotation curves of 
F568-3 for the zero, minimal, popsynth, and maximum 
disk models. We also show the constrained NFW halo 
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Fig. 13. — Isothermal and NFW halo fits for F583-1. Line types are described in Figure 5. [See the electronic edition of the Journal for 
a color version of this figure.] 



for each dark matter rotation curve. In the right panel 
of Figure 15, we plot the same data in terms of mass 
using M = V 2 R/G. To prevent a cluttered plot, we 
show only the results for the popsynth model. The solid 
line is the enclosed mass as a function of radius for the 
NFW C onstr halo. The dashed line is the enclosed mass 
of the observed dark matter. We have not smoothed the 
data, so Vdm is a bit jittery between ~ 6-8 kpc; M(r), 
however, is displayed using a smoother version of the 
data. The NFW con str halo has been forced to match 
the data at large radii; interior to this, however, there is 
less mass observed at all radii than is expected. 

In Figure 16 we plot the ratio of the observed dark 
matter mass to the constrained cusp mass as a function 
of radius for each galaxy. This ratio approaches 1 at 
large radii where the data have been forced to agree. 
Near the centers of the galaxies there is a substantial 
cusp mass excess; there is at least two times more mass 
expected in the cuspy halo than is allowed by the data. 



At all radii, the cusp mass excess becomes larger as the 
baryons become increasingly more important. 

8.2. Reconciling the Cusp Mass Excess with Noncircular 

Motions 

In our analysis so far, we have assumed that the ob- 
served velocity is the circular velocity (Vdrc = Vobsv)- If 
noncircular motions are present, the true circular velocity 
may be underestimated. This is the argument commonly 
given to explain the discrepancy between the NFW halo 
and the observations at small radii (e .g. IS waters et al.l 
I2003al Ivan den Bosch fc Swatersl l200ll ). To determine 
the true circular velocity, the noncircular motions are 
added in quadrature to the observed velocities, V 2 rc — 
V 7 2 ot + (5a 2 , usually assuming an isotropic dispersion (/? 
= 3). 

We can invert Figure 16 and determine how large 
the noncircular motions must be in order to bring the 
observed dark matter velocity into agreement with the 
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Fig. 14. — Left: Isothermal halo parameters for different assumptions about the baryons. The red pentagons are maximum disk fits 
where R c is a lower limit and po is an upper limit. The data for NGC 959 are beyond the range of this plot and are not shown. The core 
ra dius is often > 1 kpc. There is a trend towards larger R c with increasing Y t . For comparison , the solid line is the po-Rc scaling law 
of Kormcndv & Freeman (2004) and the dashed line is the po-Ro relation of Spano ct al. (2007). Right: N FW f ree halo parameters for 
different assumptions about th e baryons. No NFW fi ts could be made in the maximum disk case. The solid lines show the range of c - V200 
predicted by the cosmology of Tcgmark ct al. (2004), whereas the dashed lines are the predictions for W MAP3 (Spcrgel ct al. 2007|). For 
both cosmologies, the width of the bands is ±ler, assuming a scatter of a c = 0.14 (Bullock ct al. 2001). The data show a steeper C-V200 
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Fig. 15. — Left: Comparison of the dark matter rotation curve for different values of T* with the NFW rotation curves expected from 
cosmology for F568-3. Right: Total dark matter mass as a function of radius for T t pop- The solid line is the dark matter mass predicted 
by CDM; the dashed line is the dark matter mass allowed by the data. The two curves have been forced to meet at large radii; interior to 
this, CDM predicts more mass at all radii than is actually observed. [See the electronic edition of the Journal for a color version of this 
figure.] 



NFW cons t r velocity. In Figure 17 we plot the required 
cr(r), assuming an isotropic dispersion, for each galaxy. 
In general, most of the galaxies need about 20 km s _1 
noncircular motions at inner radii for the data to be con- 
sistent with the expectations of CDM. This is roughly 
twice as la rge as the velocity dispersion in the DensePak 
data (|K06T ). It is also important to recognize the effect 
of stars not having zero mass. Instead of just having to 
convince ourselves that ~20 km s _1 noncircular motions 
should be added to the data, we actually must be willing 



to consider noncircular motions an additional 5-10 km 
s _1 higher when Y* 7^ 0. Because we have forced the 
halo velocity to match the observed velocity at large ra- 
dius, our derived er(r) often show a steep drop to km 
s _1 from ~ 20 km s . Strictly speaking, this behav- 
ior is unphysical, and a more continuous model could be 
derived. If a contributes at large radii, however, that 
implies that the halo has even more dark matter and a 
at small radii would have to be even higher. In a forth- 
coming paper, we model the NFW halo and noncircular 
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Fig. 16. — Ratio of the observed dark matter mass to the constrained cusp mass as a function of radius. The black (solid), green (dotted), 
blue (short dash), and red/magenta (long dash/dot-long dash) lines are the zero, minimal, popsynth, and maximum disk cases, respectively. 
Near the centers of the galaxies there is a substantial cusp mass excess with at least two times more mass expected in the cuspy halo than 
is allowed by the data. [See the electronic edition of the Journal for a color version of this figure.] 



motions, and Figure 17 provides us with testable predic- 
tions. For instance, we can use the results to test whether 
the observations are consistent with NFW+a(r). 

9. CONCLUSIONS 

We have presented updated DensePak velocity fields, 
rotati on c urves and zero disk halo fits for three galaxies 
from IK06I . We have also presented the velocity fields, 
derived rotation curves and zero disk halo fits for six 
galaxies previously unobserved with DensePak. Overall, 
we find the isothermal halo to be a better description 
of both the shape of the rotation curves and the central 
densities of these nine galaxies in the limit of zero disk 
than the NFW halo. When NFW fits could be made, 
the concentrations were often beyond the range of values 
expected for ACDM. We also find that the quality of both 
the isothermal and NFW halo fits is greatly improved 
when the radial range of the data extends into the flat 



part of the rotation curve. 

For those DensePak galaxies with photometry, we have 
also presented isothermal and NFW halo fits for four as- 
sumptions of the stellar mass-to-light ratio, T*. We have 
tested the zero, minimal, popsynth, and maximum disk 
cases. We have found that the NFW halo is a poor de- 
scription of the data, and that the NFW halo and maxi- 
mum disk are mutually exclusive. There is a substantial 
cusp mass excess near the centers of the galaxies with at 
least two times more mass expected in the cuspy CDM 
halo than is allowed by the data. Most galaxies in the 
sample require ~20 km s _1 noncircular motions to rec- 
oncile the differences between observations and the NFW 
halo. Even larger noncircular motions are required when 
stars are allowed to have mass. We have also found the 
data to favor a low erg cosmology. 

10. ACKNOWLEDGEMENTS 
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Fig. 17. — The required noncircular motions, assuming an isotropic dispersion, to reconcile the difference between the N FWconstr 
velocity and the observed dark matter velocity. The black (solid), green (dotted), blue (short dash), and red/magenta (long dash/dot-long 
dash) lines are the zero, minimal, popsynth, and maximum disk cases, respectively. [See the electronic edition of the Journal for a color 
version of this figure.] 
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